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Abstract 

The gauge symmetries of a general dynamical system can be system- 
atically obtained following either a Hamiltonean or a Lagrangean ap- 
t:j- I proach. In the former case, these symmetries are generated, according 

^^ • to Dirac's conjecture, by the first class constraints. In the latter ap- 

proach such local symmetries are reflected in the existence of so called 
Tj- I gauge identities. The connection between the two becomes apparent, 

^^ ' if one works with a first order Lagrangean formulation. Our analysis 

applies to purely first class systems. We show that Dirac's conjecture 
, applies to first class constraints which are generated in a particular 

Aj' iterative way, regardless of the possible existence of bifurcations or 

r^ • multiple zeroes of these constraints. We illustrate these statements in 

terms of several examples. 
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1 Introduction 

The problem of revealing the gauge symmetries of a Lagrangean has been 
the subject of numerous investigations [1]-[13]. Local symmetries of a dy- 
namical system have been studied both within the Lagrangean as well as the 
Hamiltonian framework. 

On the Lagrangean level there exist well known algorithms [9, 10, 11, 12, 
13] for detecting the gauge symmetries of a Lagrangean. It has the merit of 
directly generating the transformation laws in configuration space, expressed 
in terms of an independent set of arbitrary functions, which leave the action 
invariant. Every one of the gauge parameters parametrizing such a local 
symmetry is directly related to a so-called "gauge identity" . The number of 
such parameters is equal to the number of independent gauge identities. 

On the Hamiltonian level the relevant action whose vanishing variation 
leads to the Hamilton equations of motion, is the so called "total action". 
The transformation laws in phase space, which leave this action invariant, 
have been conjectured by Dirac to be generated by the first class constraints 
[14]. The number of such constraints is in general larger than the number 
of gauge-identities of the Lagrangean formulation, referred to above. Hence 
the number of free functions parametrizing the symmetry of the total action 
is in general less or equal to the number of first class constraints. This is a 
well known fact [2, 5, 6]. 

In this paper we shall restrict ourselves to purely first class systems. 
Starting from an equivalent first order Lagranean formulation, whose Euler- 
Lagrange equations are just the Hamilton equations of motion, we make use 
of well known Lagrangean methods, in order to establish a direct connection 
between the gauge identities and the generators of gauge symmetries in the 
Hamiltonian formulation. In particular we show, that for the first class con- 
straints to be generators of local symmetries of the total action, these must 
be generated in a particular iterative way. It is in this sense that Dirac's 
conjecture holds. This generalizes the results of one of the authors, ref. [15], 
to an arbitrary number of primaries. We use this framework in order discuss 
several examples which have been cited in the literature as counterexamples 
to Dirac's conjecture, in order to point out some subtleties, and to confirm 
Dirac's conjecture in the sense described above. This is the main objective 
of the paper. 

The organization of the paper is as follows. In section 2 we discuss the 



iterative generation of the first class cliains terminating in gauge identities 
for an arbitrary number of primaries, following the general ideas of ref. [15]. 
We then derive from the gauge identities the transformation laws for the co- 
ordinates, momenta and Lagrange multipliers, which leave the total action 
invariant, and show that they have the form conjectured by Dirac, if the 
first class constraints are generated in a definite iterative way. We illustrate 
in section 3 the formalism for the case of an arbitrary system with two pri- 
mary constraints and one secondary constraint. Section 4 is devoted to three 
concrete examples which have been considered in the literature as counter 
examples to the Dirac conjecture: the first one discusses the possibility of 
a first class constraint becoming "ineffective" on the level of the equations 
of motion, though generating an off-shell symmetry on the level of the total 
action; it also demonstrates that the replacement of a first class constraint 
generated iteratively by an equivalent one is in clash with Dirac's conjecture. 
The second example illustrates the possibility that our algorithm generates 
the correct local symmetries of the total action, away from the constrained 
surface, eventhough the "structure functions" may be singular on that sur- 
face. The third example involves bifurcations of an iteratively generated 
constraint. We show in this case that the choice of a particular bifurcation 
will not correspond to an off-shell symmetry of the total action, while the 
strict implementation of our algorithm uncovers the full local symmetry of 
the total action. We conclude in section 5, where we once more summarize 
our general statement regarding Dirac's conjecture. 

2 Gauge identities and Dirac's conjecture 

The following analysis applies to purely first class systems, and is a exten- 
sion of ref. [15] to the case of an arbitrary number of primaries ^ . Let 
0^°-' = (o = 1,2, ■ ■ ■,A^) be the primary constraints associated with a 
seconc? order Lagrangean L{q,q), where qi,i = l,---,n are coordinates in 
configuration space, and let H{q,p) be the corresponding Hamiltonian de- 
fined on the primary surface. One readily verifies that the Euler-Lagrange 

^When refering to that reference, the reader should be aware of a number of notational 
differences. 



equations associated with the first order (total) Lagrangean 

n 

Lt(?,jo,^,P, A, A) = ^PiQi - HTiq,p,X) (1) 

with 

N 

HT{q,p,X) = Hiq,p) + J2X^<P^^\ (2) 

a=l 

reproduces the Hamilton equations of motion including the primary con- 
straints, if we regard qi,Pi and A" as coordinates in an 2n + A^ dimensional 
configuration space. We now write (1) in the form 

2n+N 

Lt= T. ^c.{Q)Qa - Ht{Q) , (3) 

a=l 



where 
and 



Q, :=(g,p,A\---,A^) (4) 



N 



Ht{Q) = H{Q,, ■ ■ ; Q2n) + E Q2n+a0i°) • (5) 

a=l 

The non- vanishing elements of a^ are given by a, = Qn+i = Pi {i = ^ ■ ■■,n). 
The 2n + N components of the Euler derivative are ^ 

dt \dQj dQ^ ^ ' 

2n+N 

= - E ^SQ/3 + ^r , (7) 

/3=1 

with 

-^0/3 = dattis — djsaa ■ (8) 



■^Our definition of tlie Euler derivative differs from that of ref. [15] by a minus sign. 
Tlris reference contains some printing errors. In particular the index "a" in Eq. (56) of 
that reference takes the values a — 2,- ■ ■ , M, and a in Eq. (59) runs over 1, 2, • • • , 2n + 1, 
and not from 1 to 7, as stated in the paragraph following (59). 



F(°) is the (2n + A^) x (2n + A^) matrix 



F(o) 



/ -10- 
10 0- 

0^ 0^ 0- 



\(F (F 0---0 J 



■ o\ 

■0 
■0 



(9) 



and 



k!"^ 



dQa ' 



(10) 



where 1 is a n x n unit martix, are A^-component Null column vectors 
(associated with the absence of A" in Lt), and 0"^ is the transpose of 0. 
The variation of the total action 



Si 



dtLT{Q,Q) 



is given by 






(11) 



(12) 



where we have dropped a boundary term. The left-zero modes of F*^") are 
given by 

v^^)[a) = (0,0,n(a)) , (13) 

where h{a) is a A^-component unit vector with the only non-vanishing com- 
ponent in the a'th place. Hence 



^/°)(a)-i?(°) = ^/°)(a)-K(°) = 0i' 



(0) 



(14) 



We thus recover on-shell the (first class) primary constraints (j)^^' = 0. 

We now adjoin the time derivative of the primaries {0*^°-*} to E^^^ and 
construct the 2n + 2N component (level one) vector E^^^ : 



^(1) 






rf.0(O) 



(15) 



where 0*^'^'' is a colum vector with components {(pi , ■ ■ ■, (p^ )■ By construction 
E^^^ vanishes on shell, i.e. for E'^'^^ = 0. The components of E^^\ which we 
label by ai, can be written in the form 



eS = -T.f2:>A + kj^^Q) 



(16) 



where 



K(i) 







and F^-*-^ is now the rectangular matrix 



/ 



F(i) 




1 

0'^ 



-1 



0'^ 



0' 
0' 



0- 



-V(/>1' 



(0) 



-V# 0- 



V-v<^S^) -WS^) o---oy 



•0 
•0 



■0 
■0 



Here 



V:=(ai,---,a„), V:=(a, 



n+l, 



.a 



'2n 



(17) 



(18) 



(19) 

We seek new constraints by looking for left zero modes of F*-"*"-*. They are A^ 
in number, and are given by 



v^^Ha) 



-V#,V0i°),O,eW(a)) , 



(20) 



where e'^'^^a) is an A^-component unit vector, with the only non- vanishing 
component in the a'th position. This leads to 



f211 



dqi dpi dqi dp, 



or 



v('\a) . i« = {0i°),/J} - EA^{#-€^} > (22) 



which by constuction vanish on shell (i.e., for E^^^ = 6). Poisson brackets will 
always be understood to be taken with respect to the canonically conjugate 
variables qi and pi. For the purpose of illustration we suppose that the second 
term on the right hand side vanishes on the surface defined by the primary 
constraints, and that we have a new constraint (p^^^ = 0, with ^ 

<Pl^^:=W^:\H}, (23) 

which is only a function of q and p. Hence from (23) and (22), 



',(1) 



.-(^)(a).i« + EA^{#,#}, (24) 



V 



or 



0« = t;«(a) . ^(^) + E X^dT{^'\b) . ^(°)) , (25) 



b,c 



where the structure functions CJ.^^ are defined by 

{#,0i°)}=E«°vr, (26) 

c 

and use has been made of (14). Note that (p^^^ is again a function of only q 
and p. 

We now repeat the process and adjoin the time derivative of the constraint 
(25) to the equations of motion to construct E^'^^: 



^(2) 






(27) 



where 0*^°-* and (p^^^ are A^ component column vectors. This leads to a matrix 
F*-^-*. As we continue with this iterative process, the number of new zero 
modes generated at each new level will in general be reduced, as "gauge 

'*If this is not the case, the algorithm stops. 



identities" are being generated along the way (see below). Hence the number 
of components of (j)^^^ will in general decrease, as the level i increases. 

The constraints c;/)^^-* with i > 1 can be iteratively constructed from the 
recursion relation, 



/,^^+i) = ^^+i)(a^ . ^(^+1) ^ Y- Y- \cr^[o«'l^(^') 



^(m) + ^^^c^[oj,^(.)^ ^>0 (28) 

£'=0 b,c 



where, 

0r^:={0i^\^} (29) 

and the sum over b in (28) runs over all constraints (pi at level £', ^ The 
coefficients CJ.^^ are structure functions defined by 

{#>i^^}=EE^rvr. (30) 

e'=o b 
The zero modes at level £ + 1 have the following generic form : 

^;(^+i)(fe) = (-Vrf),vrf\0,eW(6)), i>0. (31) 

Here e'^^^b) is a unit vector with the only non- vanishing component at the 
position of the constraint (pl in the array {(j)^'^\ <j)^^\ ■ ■ ■^(j)'^^^) appearing in 
the expression for E^^\ The iterative process in a chain labelled by "a" will 
come to a halt at level i. = Na + 1-, when 

e"^^^ = {e"\^}=EE4r'0?. (32) 

1=0 b 

Making use of (28) and (32), and setting i = Na in (28), this equation takes 
the form 

Ga := v^^'^^'Xa) ■ ^(^"+^) - E E ^If "tf ^ = , (33) 

i=0 b 

where 

i^ir' = /^ir-EA^^ir^ (34) 

c 



^Note that this set of constraints may be smaller in number than the set of primary 
eonstraints. 



Eq. (33) expresses the fact, that the 0-chains labelled by "a" ends in a "gauge 
identity" at the level A^^ + 1. 

Iteration of (28), starting with i = 0, allows us to express all the con- 
straints in terms of scalar products v^^^ ■ E^^\ Substituting the resulting 
expressions into (33), and multiplying each of the gauge identities (?„ = by 
an arbitrary function of time aa{t), the content of all the identities can be 
summarized by an equation of the form 



N Na + l 

a=l e=o 



where 






(35) 



eW 



dt 



i<t>^'' 



dt 



(36) 



Now, because of the generic structure of the eigenvectors (31) we have 
from (36) 



2ra 



v^^^^\a) ■ ^(^+^) = E {v^^'\a)E^^) + ^ , 



0, 



iV„ 



(37) 



where n is the number of coordinate degrees of freedom, and the constraints 
'Pa^ appearing on the RHS can be expressed, by iterating (28), in terms of 
scalar products v'^^^ ■ E'^^\ which in turn can be decomposed in the form 
(37). Upon making a sufficient number of "partial decompositions" udv = 
d{uv) — vdu, the identity (35) can be written in the form 



N Na 



2n+N 



a=l e=0 a=l 



,(^+l)( 



N 2n+N 



)E!^n + T.~^aT. 



v^\a 



)E^^') 



a=l a=l 



dF 
'dt 



0, (38) 



where the {e^^} and €a depend on the A^ arbitrary functions of time {^^(t)}, 
as well as on the QaS and time derivatives thereof. The expression is of the 



form ^ 

-E^i°^^Q« = f> (39) 

where 

N Na N 

SQ. = -T.T. 4'^e^H«) - E ~^y'\^) ■ (40) 

a=l e=0 a=\ 

For infinitessimal e^^-* the time integral of the LHS of (39) is just the variation 
of the total action (12). Hence we conclude that the transformations (40) 
leave the total action (11) invariant. But because of the generic structure of 
the eigenvectors (31) we have from (40), 

^?. = ^Q^ = Y.Y.^T-^ = 1:Y. 4^H^. m ,^ = l,■■■,n. (41) 

a=l e=o ^P'' a=l e=o 

5p. = 5Q^^. = - E E 4'^^ = E E 4'Hp. <^i^)} , (42) 



=1 e=o "^y* a=i e=o 



while 



.^A-^ = dQ2n+a = -ea • (43) 

Note that according to (29) the constraint cp^^'^^' has been generated from cp^^' 
via the Poisson bracket of (j)^^^ with H. Hence, given the primary constraints, 
the form of the secondary constraints are fixed, and one should not modify 
their structure, replacing them by an equivalent set of constraints. ^ We 
have thus proven a proposition made in [16] regarding the form of the first 
class constraints that are generators of local symmetries. 

Finally we would like to emphasize, that we have implicitely asumed 
throughout that the structure functions appearing in equations (30) and 
(32) are finite on the constraint surface. If this is not the case, we expect 
that our algorithm will still generate the correct local off-shell symmetries of 
the total action, if the gauge identities are non-singular on the constrained 
surface, while they do not implement symmetries on the level of the Hamilton 
equations of motion. This will be demonstrated in our example 2 in the 
section 4. In the following section we first illustrate our formalism for the 
general case of two primary constraints and one secondary constraint. 

^Thc minus sign has been introduced to cast the transformation laws in a standard 
form, when written in terms of Poisson brackets. 

^Thus, for example, one should not replace a constraint p^ = by p = 0. 



3 General system with two primaries and one 
secondary constraint 

In this section we apply the formahsm to the case where the system exhibits 



two primary constraints 



, a 



1,2, and one secondary constraint 



lW 



which in Dirac's language is generated from the presistency in time of 02 • 
(This is illustrated by the table below.) The specific form of the Hamiltonian 
is irrelevant. The constraints are assumed to be first class. 



^(0) 



4°) 



^ 



(44) 



The total Lagrangean L^ is given by (1) with N = 2. Proceeding as in 
section 2, the primary constraints can be written in the form (14), where 
y(o)(l) = (0,0,1,0) and vW(2) = (0,0,0,1). 

Next we construct the vector (15), i.e. E^'^^ = {E'^^\(j)\ ,02 )i ^"^^ the 
zero modes of the corresponding matrix F*^^\ which are given by (20), i.e., 
v^^H) = (-V0S°\V0S°\O,O,1,O), v('\2) = (-V0f ,V#,0,0,0,1). Since 
only 02 is assumed to lead to a new constraint, we are immediately left with 
one gauge identity at level 1, generated from 0^ : 



or 



Gi=.-(^)(i)-^«-E^r0r 

b=l 



Gi = ^(^)(i) ■ ^(^) - y: ^r (^'\b) ■ ^(0) 



(45) 



(46) 



(0) 



which is nothing but (33) with A'^i = 0. On the other hand 02 gives rise, by 



(1) 



assumption, to a new constraint 02 at level 1, which is given by 



^(1) = ,7(1) 



(2)-i«+ ^ A^cSr(^;(°)(fe)-i(°)) 



(47) 



We are therefore led to construct 



m 



i,C=l 



A^(o) 
A^(o) 



(48) 



10 



as well as the corresponding rectangular matrix F^"^^ and its left zero modes, 
of which only the contraction of 

t;(2)(2) = (-V#,V0(^\o,O,O,O,l) (49) 

with E^"^^ leads to a new equation, which is necessarily a gauge identity, since 
we have assumed that the system only possesses one secondary constraint 
02 • The gauge identity at level 2 has the form 

G, = v('H2) . E(^) - ± 41°' ^r - Kg''<P'2' - , (50) 

b=l 

or 



G2 = v^'K2) ■ E^'^ - ^4l°](t;(o)(6)-^(o)) (51) 

6=1 



6,c=l 



Multiplying the gauge identities (46) and (51) by the arbitrary functions ai{t) 
and a2{t), respectively, and taking their sum, the resulting expression can be 
written in the form (35). Upon making a sufficient number of "partial differ- 
ential decompositions", one finds, after some algebra, that the information 
encoded in the gauge identities can be written in the form (38), where 

(0) 

ey = as + a2-f^22 

6x^ = -d2EA^«°'-«24i°'-«i4r 

c 
- «24?EA'^Sa°'+«252a-ai^la + ^(«242'V2- ^^2) 

are the only non-vanishing parameters. The corresponding transformation 
laws for the coordinates qi and momenta pi are given by (41) and (42). 

Having obtained e^^ , £2 , €2 and SX"" expressed in terms of ai{t) and 
a2{t), one now verifies that the e^^-* are solutions to the well known recursion 



11 



relations following from the requirement that the transformation laws (41) 
and (42) be symmetries of the total action[2, 5, 6]: 

4^^ + EE4'^^r = o, ^>i (53) 

£' b 

and 

«"=4°'+i:i:er'Ar^ (54) 

e' b 
In our case these equations reduce to 

^ (1) 2 

^ + E4°^<^+4^^<^=0, (55) 

"■^ b=i 

and 

^A« = ei°) + Eer4r+4^^4i^ (56) 

b 

where we have made use of the fact that, because of the way in which the 
constraints have been generated, K^i = (a = 1,2), K22 = 1, Ku = 0, 
Cabc = 0, and K22 = —Y.c^'^dc22- One then verifies that the above 
equations are indeed satisfied. 

4 Counterexamples to Dirac's conjecture? 

In the literature it has been stated that Dirac's conjecture is not always 
correct [17, 18, 19]. In the following we present three models which have 
served as examples for this in the literature and show, that there is no clash 
with Dirac's conjecture if the first class constraints are understood to be 
generated in the definite iterative sense discussed in section 2. 

Example 1 

An example considered in [18] is given by the Lagrangean 

L=^e'^m (57) 

In the following we analyze this system in detail, and point out some sub- 
tleties leading to an apparent clash with Dirac's conjecture. 

12 



The Lagrange equations of motion can be summarized by a single equa- 
tion 

qi = 0. (58) 

Hence g2 is an arbitrary function. 

Eq. (58) does not possess a local symmetry. On the Hamiltonian level 
the system nevertheless exhibits two first class constraints, which, as we 
now show induce transformations which are off-shell symmetries of the total 
action J dt L^, with L^ defined in defined in (1). However only one of the 
constraints generates a symmetry of the Hamilton equations of motion. 

From (57) we obtain for the (only) primary constraint 

0(0) = p2 = . (59) 

The canonical Hamiltonian evaluated on the primary surface is given by 

if(g,p) = ie->^ (60) 

and the first order total Lagrangean reads 

Lriq, P, A; g, p,X) =J2 PiQi - H{q, p) - A(/)^°) . (61) 

i 

Considered as a function of the "coordinates" qi, pi and A and their time 
derivatives, it has the form (3), where 2?7, + A^ = 5, Qa = {q,p,\), and 
c^i = Pi (^ = 152) for the no n- vanishing components of a^. 

The Euler derivatives are given by (6), where F^^ = d^ap — dpUa are 
the elements of a 5 x 5 matrix with non-vanishing components F13 = F24 = 
— F31 = — F42 = —1. The equations of motion are given by £"^0) = q^ ^j^d 
yield the Hamilton equations of motion 





Qi - 


- e"' 


''Pi = 


(12 - 


- A = 


= 


Pi = 


= 

1 




P2 - 


- r 


--P? 



0, (62) 

as well as the constraint (59). We now proceed with the construction of the 
constraints and gauge identities as described in section 2. Since we have only 
one primary constraint, the formalism simplifies considerably. 

13 



The matrix F^^ has one left zero-mode v^^^ = (0,0,0,0,1), whose contra- 

tion with E^^^ just reproduces on shell the primary constraint: ■u'-'^-' ■ E'^'^^ = 
0(0) _ Proceeding in the manner described in section 2, we construct E'^^^ and 
corresponding left eigenvector of F^^\ v'^^^ = (0, —1, 0, 0, 0, 1), leading to the 
secondary constraint -y*-^-* ■ E^^^ = (f)^^\ where 



id) 



dt 



m 



dV 



7(0) . 



^(0)) 



and 



0(1) 



1 



e~«2p^ 



(63) 



(64) 



Note that by construction cf)^^^ = 0, on shell. Since {(j)^'^\ cp^^^} = 4'^^\ 
and 0(1-' form a first class system. The algorithm is found to stop at level 
"2" where the following gauge identity is generated: 



with 



and 



^2) 



^(2) 



^2) . ^(2) ^ ;^^1) . ^(1) _ Q 
1 



-e'>i,0,0, 



'92^2 



Pi, 0,0,1) 



/ F(0) \ 



The gauge identity can be reduced to the form 



dt 

d<t>W 

dt 



( 



ew 



\ 



v^-' ■ F(o)) 



dt 



(65) 
(66) 

(67) 



-f i?f + |(^^«i?f ) + ^(-f 4°)) + A[(.;i^)i?(°)) + |(.;^ (68) 

where a summation over a = 1, ■ ■ ■, 5 is understood. Multiplying this expres- 
sion by an arbitrary function of time e{t), this identity becomes of the form 
(39), with 

d 



SQa = e<^ + (Ae - e)vl'> + e 



(1) 



dt 



(A. 



,(0) 



(69) 



In terms of the Hamiltonian variables, g^, Pi and Aj, (69) implies the following 
transformation laws 



6qi 



—ee '^^pi 



6(^H?i,r^}, 



1)~ 
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Spi = 0, 

5X = ^■-|(Ae)=e(^ (70) 

where 

e(0) = e-Ae, e« = -e . (71) 

One readily verifies, that the e'-^^ satisfy recursion relations (53). 

We have thus verified that the transformations correspond to a symmetry 
of the action (11). This symmetry is realized off shell, and requires the full 
set of transformation laws (70). On the other hand, the first class constraint 
^(1) := p-^ = 0, although equivalent to the secondary constraint (f)^^' = 0, is 
not a generator of a local symmetry of the total action. Indeed, as observed 
in [18], pi induces translations in qi, which is not a symmetry of the total 
Lagrangean (61). Note further that on the constrained surface, the variations 
of dqi and dp2 vanish, and hence do not generate a symmetry on the level 
of the Hamilton equations of motion. Thus 0*^^-* becomes and ineffective 
generator on this level. The remaining symmetry is just the statement 6q2 = 
dX, which is consistent with (62). It has no analogue on the level of the 
Lagrange equations of motion following from (57). 

Example 2: constraint of the form (p = f'^ = 0. 

The following example taken from ref. [19] illustrates our comment made 
at the end of section 2 regarding the case where the structure functions (34) 
are singular on the constrained surface, while the gauge identitites are well 
defined on that surface. 

Consider the Lagrangean, 

L = ^e + uf\q), k>l, (72) 

which has been classified as being of "type II" in ref. [19]. ® The canonical 
Hamiltonian evaluated on the primary surface Pu = is given hj H = ^p"^ — 



^In order to simplify the discussion we have restricted ourselves to two degrees of 
freedom, q and u. 
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uf^{q), and correspondingly the first order Lagrangean reads 

\ 
MQ, Q)=PQ + PuU - -p^ + uf^{q) - Xpu , (73) 

where Q = {q,u,p,Pu, X). The Euler-Lagrange equations derived from (73) 
are just the Hamilton equations of motion (including the primary constraint 
Pu = 0). Following the above Lagrangean algorithm, we obtain at levels 
£ = 0,1,2, 



V 



(0) _ m n n n 1^ J^) 



V 



(0,0,0,0,1), 0W = ^/o).^W=p„, (74) 

« = (0,-1,0,0,0,1), 0« = ^/i) . E« = /^(g) , (75) 

(2) = (0,0,r-Y,0,0,0,l), v^'^-E('^ = kp{\nfycf>('\ (76) 



where the prime denotes the derivative with respect to q. Note that for k > 1 
the right hand side of (76) is well defined on the constrained surface (p^^' = 0, 
and in fact vanishes there. Hence at level £ = 2 we have generated a gauge 
identity, which upon making use of 0'-^^ = ■u'-^-' ■ E'^^^ takes the form. 

In the present example the structure functions in (30) vanish identically, 
and the "coordinate" dependent factor kp{lnf)' in (77) is nothing but the 
structure function K in (34). This function is singular on the constrained 
surface. Nevertheless the gauge identity (77) "generates" Oj^-s/ie// the correct 
local symmetry of the total action, as we demonstrate below. 

In the present example E'^'^^ and E^"^^ have again the form of (63) and 
(67), respectively. Making use of these relations, and multiplying the gauge 
identity (77) by e, one finds that the resulting identity can be written in the 
form (39) with dQa the infinitesimal transformations, 

(78) 



(79) 



Sq = 

6u = 

6p = 

Spu = 


0, 

e + ekp{lnfy, 

ekf-'f , 

0, 


6X = 


e+|(6Mln/)')- 
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These transformations are only defined off-shell. Computing the correspond- 
ing variation 5Lt one finds that it is given by a total time derivative ^ : 

5Lt = j^{ef\q)) . (81) 

In accordance with our expectations, one verifies that 5Qa can be written in 

the form 

1 

<5Q« = Ee^{Q"'</'^'^}' « = l,---,4, (82) 

where 

eo = e + ekp{\nf)', (83) 

ei = -e. (84) 

For the Lagrange multiplier one has (5A = J^eo- It is then an easy matter to 
verify that e^ and 5\ indeed satisfy equations (53) and (56). 

The above infinitesimal transformations (78) represent symmetry trans- 
formations of the total action, away from the constraint surface f^{q) = 0. 
They are not defined on the surface (p^^^ = 0. Correspondingly they do not 
represent symmetries of the Hamilton equations of motion. 

Summarizing, the above example shows that if we generate the constraints 
in a strict iterative way following the Lagrangean algorithm, then the first 
class constraints do generate the local off shell symmetries of the total action 
in the form conjectured by Dirac, if the gauge identities are finite expressions 
on the constraint surface, although the structure functions may be singular 
there. We have, however, no proof that this is generally true. 

Example 3: bifurcations 

In the following we now consider an example of a system exhibiting bifur- 
cations of constraints, such as considered in ref. [20], and show that also in 
this case our algorithm generates correctly the gauge symmetries of the total 



^For the second order Lagrangean the corresponding gauge identities are found to imply 
the infinitesimal symmetry transformation, 

dq = 0, 6u = e + ekp{ln /)' . (80) 

One again verifies in this case that SL — ■§^{e.f''{q)) . 
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action, whereas the first class constraints corresponding to the choice of a 
particular branch of these bifurcations are not part of our iterative algorithm, 
and do not generate gauge symmetries of the total action. 

Consider the so called "Christ-Lee model" [21] discussed in ref. [20], and 
defined by the Lagrangean 

L = -ql + -ql - qMq2 - q2qi) - V{ql + ql) + -ql{ql + ql) , (85) 

with V{x) some potential. The equivalent first order Lagrangean in our 
approach reads, 

Lt = ^Piqi - ^pI - 2^2 - q3{qiP2 - q2Pi) -V -\p^. (86) 

i 

The corresponding Euler-Lagrange equations take the form (6) where Qa '■ 
{qi,q2,q3,Pi,P2,P3,>^) and the E^^^'s are given by 



,(0) _ 

1 


= pi+P2q3 + diV, 


Ei'^=P2-piq3 + d2V, 


i{0) 
'3 


= P3 + {qiP2 - q2Pi) 


, Ef^ = -qi +pi- q2q3 


,(0) 
'5 


= -q2+P2 + qiq3, 


^6 - -q3 + A , -C/7 - 



P3- 

Going through the procedure of section 2 we arrive at a level zero and level 
one constraint, 

(f)(0) ■- yio) . e(^) = p^ = ^ 0(1) :={;(i).^(i) = g2Pi-giP2 = 0, (87) 

where v^^^ = (0, 0, 0, 0, 0, 0, 1) and v^^^ = (0, 0, -1,0, 0, 0, 0, 1), and E^^^ has 
the form (63). The iterative process stops at level two, with the gauge identity 
^2). ^(2) ^Q^ ^(2) has the form (67) with t;(2) = (-g2, ?i, 0, -p2, Pi, 0, 0, 0, 1) 
the level 2 zero mode. The simple form of this identity refiects the fact that 
the structure functions K^f^ appearing in (34) vanish for the case in question. 
Explicitely one has 

,K2) . ^(2) ^ ^(2)^(0) + 1^(1)^(0) + ^.(0)eW , 0. (88) 

This equation, when multiplied by e(t), can be written in the form (39) with 
F = e(j)^^^ — e(p^^\ The {6Qa} can be compactly written as 

1 
5qa = Y.^dqa,<P^'^}, S\ = e, (89) 



e=o 



U 



where eo = e, ei = — e. 

The constraints (p^^^ = and (p'^^^ = are first class and are thus found to 
generate an off-shell gauge symmetry of the total action, in agreement with 
Dirac's conjecture. 

Suppose now that we had chosen instead of (p^^^ = any one of the 
bifurcations: i) Qi = q2 = 0, ii) Pi = P2 = 0, iii) Qi = Pi = 0, and iv) 
12 = P2 = 0. Consider for instance the case i) ^^ . Then we continue 
at level one with the constraints, (pl := qi = (p2 := q2 = 0- The 
requirement of persistence in time then leads to two new constraints at level 

—(2) ^(2) —(2) —(2) 

two: (pi = Pl = , (p2 ■= P2 = 0. Since (pi ~ and (p2 ~ there are 

no further constraints. We thus see that the system of constraints has been 

reduced to a mixed system with two second class constraints, and just one 

first class constraint (p'^'^^ '■= Ps = 0. Correspondingly we are led to consider 

G = ep3 as a potential generator for a local symmetry. One readily verifies 

that all variations vanish except for dq^ and 6X. Hence 

SLt = (e - S\)p3 - e{qiP2 - q2Pi) (90) 

To have asymmetry we must demand that 5\ = e and qiP2 — q2Pi = 0, so 
that Lt only exhibits a local invariance on the restricted surface defined by 
the first class constraint generated at level 1 by the algorithm of section 2. 
Hence G does not generate an off-shell symmetry of the total action. The 
Hamilton equations of motion, however, exhibit a symmetry Q's -^ Q's + e, 
A -> A + e. 



5 Conclusion 

In this paper we have first generalized the work of ref. [15] to arbitrary first 
class Hamiltonian systems. Using purely Lagrangean methods we have de- 
rived the local symmetries of the total action, and have discussed in which 
sense the Dirac conjecture holds. More precisely, the Dirac conjecture was 
found to hold for first class constraints that have been generated in a definite, 
strict iterative way by the Lagrangean algorithm described in section 2. This 
provides a proof of the prescription given in ref. [16] regarding the form of the 

^°We are leaving herewith our algorithm, since with the choice of a particular branch 
the constraint can no longer be written in the form "tr^' ■ E^^' . 
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first class secondary constraints generating a local symmetry. In particular, 
we have quite generally shown that the infinitesimal gauge transformations 
which leave the total action invariant have the form proposed by Dirac, i.e., 
^Qi = T,A^^{QiAA}, Spi = T,A(^^{PiAA}, where the {(/)a}'s are all the first 
class constraints generated by the algorithm refered to above. On the other 
hand, the replacement of constraints generated iteratively by an equivalent 
set of constraints is, in general, not allowed, and in fact may obliterate the 
full symmetry of the total action, as is demonstrated by example 1 in section 
4. Our example 2 in that section illustrated that our algorithm generates cor- 
rectly the symmetries of the total action away from the constrained surface, 
although the structure functions may be singular on that surface, provided 
the gauge identities are well defined everywhere. However, for the Hamilton 
equations of motion to exhibit this symmetry, the structure functions must 
themselves be well defined on the constrained surface. Finally, example 3 in 
section 4 illustrated that in the case of bifurcations of constraints, the choice 
of a particular branch may change the structure of the first class constraints, 
such as to be in confiict with Dirac's conjecture. We emphasize that such a 
procedure is excluded in our algorithm. 
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